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INVERSE L TO I TRANSFORMATION 


SUMMARY : It is shown that the second adiabatic invariant I may 

be defined as a function of the magnitude of the field B and 
the shell parameter L. A method is developed to obtain 
solutions for I in terms of arbitrary B and L. Results 


are tested and discussed. 



INTRODUCTION 


The wide acceptance of Mcllwain's magnetic coordinates B and L 
(Mellwain 1961) and their almost exclusive use in the mapping of geo- 
magnetically trapped particles, has resulted in the gradual disuse of the 
second adiabatic or longitudinal invariant I. Many studies pertaining 
to the motion or the spatial distribution of trapped particles, employ 
predominantly the parameters B and L in developing their theories 
(Galbraith 1965), analysing phenomena (Hess 1964 ), and presenting ex- 
perimental data (Mellwain 1963 ) • 

Already a great amount of data have been given exclusively in L 
without any reference to I. This trend will probably prevail even more, 
as new and faster methods of obtaining L are being devised and used 
(Stem 1965 )j making its calculation entirely independent of I. 

When considering however topics like the distribution, the shell 
configuration or the drift of these particles, it might be desirable, 
if not advantageous, to use I in addition to B and L or instead of L. 

In most such instances the integral invariant will not be readily available 
and will have to be obtained by same computational method. 

The conventional computation of I may not be of help^be cause it can 
only be performed, if the appropriate geographic coordinates are known, in which 
case one must resort to the relatively slow process of integrating numeri- 
cally along a field line an expression equivalent to that defining I: 

1 * 

where B(s) is the field strength along a line of force, B and B* are the 

m m 

magnitudes of the field at the conjugate mirror points, and ds is the 
differential path length on a field line from B to B*. Also it is not unlikely, 

in m 

that sometimes a line integration to determine L has already been performed and 
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we just wish to save another integration when deriving 'I'. 

This paper will present a simple way, by which 'I' may be obtained 
from the known quantities B and L only, without the need of integration. 



METHOD 


Let us define the integral invariant as a function of the magnetic 
parameters B and L: 

I = f(B, L) (2) 

and let us stipulate that an acceptable solution to this function should 
be positive, real and single valued. 

Recalling from Mcllwain (1961) that: 



where : 


F (ir) - 1 + ex * (ps/) 


and: 



m = 0, 1 ... 9 (4) 


(5) 


♦Equation (4) is an approximation to a high order of accuracy, where the 

coefficients a have been evaluated for different ranges of 
m 

a (Table l). 
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we require that a solution to (2) must he derivable from (3). 

! Writing out equation (3) and rearranging, we get: 

ln (tt ' 0 = sy* m s e < 6 > 

By defining p as a function of o = g(l,B) we have reduced the problem 
to the level of plain matrix manipulation, where: 


m 

Pi “ 


m = 0, 1 . . .n 

1-1, — , n + 1 (7) 


n = 9 

will give us a system of i = n + 1 equations, for which the values of 
a may be autonomously generated, that is, arbitrarily assigned in such a 
way, as to be evenly distributed over the range of the given set of 
coefficients a^. Since equation (7) is monotonic within the fixed range, 
the g_Js are nonsingular. Hence it can be shown (appendix) that for all 
pairs of-, (p , ^ ) of a domain, some new set of coefficients b^ exists, 
that will make the inverse relation hold true: 

m = 0, 1 ... n 

Oj = i = 1> — , n + 1 (8) 

n = 9 

Combining ( 5 ) and (8) and writing the equation in a farm corresponding 
to (3) we have: 
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We shall call this expression the "inverse Dipole Function" and define 
it as: 



where : 

F ' 1 (ir) “ ex P m = 0, 1 - 9 (11) 

and: 

e - Crr - 0 (12) 

Expression (2) is now sufficiently determined from (10) by: 



To obtain the new coefficients b^, the corresponding p's were computed 

through equation (7) for a selected set of (n + l) a 's, belonging to 

a given range of a^. Using these values of a and p, the resulting 

linear system of equation (8) was solved for b^. Table 2 gives the 

values of b calculated by this method for the indicated domains of a . 

m m 

CONCLUSIONS AND RESULTS 

We have shown that a solution to I = f(B, L ) , equation (2), exists 
in the form of the Inverse Dipole Function, which does satisfy the 
imposed requirement of derivability . When conforming with the expressed 
stipulations, the second adiabatic invariant may unambiguously be 
obtained for any arbituary B and L with great accuracy and speed and 
without integration, from: 
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where : 

F_1 (ir)" eXp fovP*) m = 0, 1 9 (11) 

and: 

b-^Ctt -1 ) (12) 

In similarity to (4), Equation (ll) is also an approximation, 
but with a higher degree of accuracy. Tests and re computations of a 
and L from the calculated 'I* have shown the average error to be less 
than 10 6 with a maximum around 10 4 . 

Sample computations are given in Table 4. They were handled by a 
program written in Fortran IV for the IBM 7094, Moonlight System. 

Actual running time per 'I' computation accounted to a very small fraction 
of a second in all domains. 
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APPENDIX (I) 




Given a polynomial of the nth degree in a that is monotonic in the interval 
L a *, for which the coefficients exist: 

n 

(A) 3 = g( a ) =2 a m a, m n = 9 

m=o 


we select (n+l) points a^, ... } such that: 

loj c^a*, an 

and calculate from (a) the pairs: 


i = 1, n+l 


3q “ s(a-|_) )j (ccq* 3 2 ~ ) > ••• 3 n+1 “ s ^ a n+ i^ ^ 

This gives (n+l) points lying on g(a). 

Further, from the monotonicity of g, it can he expected that, if the q^s 
are distinct, so are the g_Js. 

Now let us write (B) as: 

(c) (g-^, ($ 2 * *'* ^n+1* ^+1^ 

and regard '(c) as giving us (n+l) points on the function: 


(D) a = g" 1 (g) 

Accordingly, this is a unique polynomial of nth degree passing through these 
points, i.e.: 

(E) a = s g m = g -1 (g) n = 9 

m=0 


Now if we change our selection of the > we could expect 

very nearly the same polynomial in (E) again, provided the 's were well 
dispersed over the interval [ a *, g 1 *] . 
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Of course, if one selected the q^s all clustered together, at one end of the 
interval L a *, a **J for example, the g^s would be clustered together also, 
and some difficulty might be experienced in solving for the r's, which in 
that case might be different. But with a reasonable selection of the foJs], 
the values for fp ] will be approximately the same, their difference 
remaining insignificantly small. 
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Me II' wain's Polynomial Expansion Coefficients for the Dipole Function Used in the 
Calculation of L (Source: Mcll wain's INVAR Subroutine) 


? 


SV 

<8 

O 

3 0 

0 o 

• • 

A i — I 

1 


oooooooo 


*ja 

o 

& 

& 


CVJ 


cvj to in 

i i 


sv 

VO c— 
00 H 

g| 

O H 


A 

CVI 

8 

S 

CVJ 


o 

I 


O 

o 

o 
1 

o 

1 

g 

>< 

4 

1 — 1 

A 

-4- 

p * 


CVJ 

CVI 

8V 

o 

31 

8 

P 

p. 

H 

S 

VO 

CU 

to 

CVJ 

6 

i 

o 

o 

f 

o 


o o o 


H CVJ rr\ 4 4- VO ^ ^ 

i i i i i i i i 

OOOOOOOO 


d d d 

O ON S S 

in cvi ^ 




H CVJ CVj' 


oi m >■ 

• • • 

o o o 

i 


0 

1 



0 o 

1 


I 

o 


*A 

V? 

« 

v 

*8 


o CO 

3 % 

VQ CVJ 
ao A 
CVJ A 
CVJ A 
VO -4* 



H CVJ 4* 

I I I 

o o o 

s ^ ll 



** 8 


in m qj cu 
H H u 5 rA 



oooooooo 


VO A A 
tA lA fc- 
VO CVJ ^ 
Vp fA IA 
O CVJ A 

3 8 8 


*8 

$ 


8 J 8 

<— I A 
CVJ A O 
VQ A 
O A 
O A 
• • 

K\ O 


O O O O O O O 


*° cp sP «r «f «r ^ «r <s° 



NEW POLYNOMIAL EXPANSION COEFFICIENTS FOR THE INVERSE DIPOLE FUNCTION USED IN THE CALCULATION 
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